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j> . Abstract 

CN . We study weak-coupling perturbation expansions for the ground-state energy of 

the Hamiltonian with the generalized spiked harmonic oscillator potential V(x) = 
Bx 2 + -4 + , and also for the bottoms of the angular momentum subspaces labelled 
by I = 0,1,..., in N -dimensions corresponding to the spiked harmonic oscillator 
potential V(x) = x 2 + where a is a real positive parameter. A method of Znojil 
Q_i' is then applied to obtain closed form expressions for the sums of some infinite series 

^ <— | whose terms involve ratios and products of gamma functions. 
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I. Introduction 

The spiked harmonic oscillator Hamiltonian defined by 

d 2 A 

H = -—+x 2 + —,a<5/2,xe[0,oo) (1.1) 

where the positive parameter A measures the strength of the singular term, has been 
the subject of intensive study [1-3]. Aguilera-Navarro and Guardiola [1], employed a 
resummation technique to obtain a weak coupling perturbation expansions for the 
ground state energy of the Hamiltonian (1.1), using standard perturbation theory up 
to second order. The Hamiltonian (1.1) is first written H = H Q + XV; then, using the 
odd-parity solutions of the one-dimensional harmonic oscillator ip n (x) = \n> satisfying 
the Dirichlet boundary condition V(0) = 0, and the unperturbed energies E n = 3 + 
An, n = 0, 1, 2, . . . , they found that the weak-coupling expansion for the ground state 
of H to the second order in V becomes 

E = E Q + \<0\x- a \0> + A 2 ]T l<0 jT a| " >|2 +-..,«< 5/2, (1.2) 

n>l 

where <0|x~"|n> given by 

.„,..-„,.. (-»)■ r(Wt+n) n = 0,l (1.3) 

V(2n + 1)! r(|)r(f) 

In their pertubation treatment, Aguilera-Navarro and Guardiola construct a function 
F , expressed in terms of a particular form of the Gauss hypergeometric series; the 
function F was found to be useful for obtaining analytic approximations for the 
ground state energy of the for the 'non-supersingular' cases: a = \ , a = 1 , and a = § . 
They also provided a weak coupling expression valid for the case a = 2. Estevez- 
Breton et al [2] derived an exact analytical result, by using the function F, valid for 
the special case a = 2 . Later Znojil [3] derived the same result by an elegant and 
economical method. It is perhaps worth noting here that in all these works, although 
the conclusions and the results were correct, there remained an error in the F formula, 
Eq.(14) in Ref.[l], deduced by Aguilera-Navarro and Guardiola, which should read 



F = 



^if-i.f-i^ij-i-^-n- 



In Sec. II of the present paper we generalize the weak coupling expansion (1.2) to 
study the generalized spiked harmonic oscillator Hamiltonian 

H = Ho + XV = -^ 1 +Bx 2 + 4 + 4' «<5/2, x£ [0,oo). (1.4) 

where A and a are positive parameters. We show that the weak coupling expansion, 
in this case, is given by 

E - 2VB 7 + BtHff A - 4 F 3 (l,l,f + l,f + l; 7 + l,2.2;l) + ..., 
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where 7 = 1 + ±V1 + 4^4 . This expression is valid for all values of a < 7 + 1, including 
a = 2. This formula allows us to obtain perturbation expansions for (1.1) and (1.4) 
valid for the bottoms of the angular momentum subspaces labeled by Z = 0, 1, ... in 
N -dimensions. In Sec. Ill, we adopt the constructive approach of Aguilera-Navarro 
and Guardiola to generalized the function F which we then use to obtain a sum for 
the infinite series 



(SL 

E 



a \2 
2 >n 



^ An{n+ 1) (7)„ n\ 

n>l 

In Sec. IV, by employing Znojil's technique [3], our generalization turns out to be 
useful for obtaining closed-form sums for other interesting infinite series whose terms 
involve ratios and products of gamma functions. 

II. Weak coupling expansions 

Recently, we have obtained expressions [4] for the singular-potential integrals 
<m\x^ a \n> of the Hamiltonian (1.4) using the Gol'dman and Krivchenkov eigenfunc- 
tions [5] 

' 1> n (x) = \n> = C n x^ 1 +^ T + m ^e-i^ x \F 1 (-n, 1 + ^VT+lA; VBx 2 ); 

2B ! + !^l+4A r(n + j + lyT+43) (2J) 
C A — - 4— - n — 1 9 

n![r(l + Iv / r+43)P ' n ~ U ' i '^"" 

where 1F1 is the confluent hypergeometric function [9] 

i F i( a 'M = E^' (a)k=a(a + l)...(a + k-l) = ^±^. (2.2) 

and the exact eigenenergies E n = \/~B(An + 2 + \/l + AA), n = 0, 1, 2, . . . , for the singular 
Hamiltonian 

d 2 A 

Ho = --^+Bx 2 + —, B>0,A>0, xe [0,oo). (2.3) 
Hall et al found [4], for a < 27, that the matrix elements <m\x~ a \n> are given by 

<m\x- a \n> = {-\) n + m B a / A 



T(7 + m) 



n\m\ T(7 + n) 

(2.4) 



y f 1lfc /m\ r(fc + 7-f)r(f-fc + n) 



X 

fe=0 



in which each element has a factor which is a polynomial of degree m + n in a . The 
relevant matrix elements <0\x~ a \n> are given [4] by 



<0\x-<*\n> ( irn°'*J F(7) r ^-D r (f+") „ _ l 2 (2 5) 

<o\x \n>-(i)B W n!r(7 + n) r( 7 )r(f) >»-°>M,--., i 2 - 5 ) 
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Thus, writing (1.4) as H = H + XV and using the standard perturbation theory to 
the second order, we find that the weak coupling expansion (1.2) now reads 

E = 2VB 1 + Bf T{ Z~p \-\ 2 B^ T2{ T l;P V - { p} , +..., «< 7 +l (2.6) 
T( 7 ) T 2 ( 7 ) g4n(7)„n! 

where 7 = 1 + \^J\ + 4A . We observe that the ratio of the n th and (n + 1) th terms of 
the sum in the coefficient of A 2 in (2.6) is 

<0\x- a \n> 2 /(E n -E o ) 7 + 2-a 1 

= H \-o( — ), as n^oo, 



<0\x- a \n+l> 2 /(E n+1 - E a ) n n 2 ' 

so that, by Raabe's test [6], this sum is convergent for a < 7 + 1 . The expressions 
(1.2) and (2.6) are accurate for A small compared to unity. It is interesting that the 
sum of the infinite series in the A 2 coefficient can be computed exactly for arbitrary 
values of a and 7 satisfying a < 7+1 . We express the sum in terms of the generalized 
hypergeometric functions p F q defined [7] by 

rp / 00 a \ {0t\)n{0t2)n---{0tp) n Z n 

p F q {a 1 ,a 2 ,...,a p ;0 1 ,fo,..-,0 q ;z) = 2^ ta\ ia\ Tr~\ — ' ( 2 - 7 > 

\Pl)n{P2) ■ ■ ■ \Pq)n n\ 



n=0 



Indeed, the sum in the A 2 coefficient implies 

V _ v (n-l)!(f) 2 

4n (7)n n! 4 (7)n (n!) 2 

_1 (n!) 2 (|) 2 +1 



4 ^o ((^ + 1 ) ! ) 2 (7)n+i n\ 

16 ^h ( 2 )n (7 + 1)™ n! 
a 2 ^ , _ a a 



(2.6 



= 16^ 4^3(1,1, 2 + b 2 + 1:7+1,2, 2; 1) 

where we have used the Pochhammer identities (z) n+1 = z(z + 1)„ and n! = (1)„ . 
This expression can easily computed for arbitrary values of a < 7 + 1 by the use, for 
example, of Mathematica. The weak coupling expansion (2.6) now reads 

E ^ +B ' T -^T^ B tf <*(U.f +i;7+i,a.*i) + .- (»») 

The results of Aguilera-Navarro and Guardiola for the special case B = 1 , A = or 
7 = 3/2 , and for the values of a = |,Q! = l,a=|, and a = 2 follow immediately 
without the necessity of special treatment for the case of a = 2 as suggested before 
by many workers in the field [1-3]. The expression (2.9) can be further generalized 
to apply to the ground state eigenenergy at the bottom of each angular momentum 
subspace labeled by I = 0, 1, 2, . . . in N -dimensions: we just need [8] to replace A with 
A -> A + (I + l/2(N - 1))(Z + 1/2(7V - 3)). For the spiked harmonic oscillator potential 
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(1.1), we set A = or we replace 7 with I + ^ to obtain a weak-coupling expansion 
valid for the bottoms of the angular-momentum subspaces in N -dimensions. 

III. The F function 

Although, our results in section (II) cover all the cases for a < 5/2 for the Hamil- 
tonians (1.1) and (1.4) the constructive approach of Aguilera-Navarro and Guardiola 
allows us to obtain more sums of infinite series involving gamma functions. We first 
generalized the function F as introduced in Ref. [1] (we also point out the error in 
the F formula there). If we denote the sum in the A 2 coefficient, of Eq.(2.6), by 



g = 2 

4-; 4n ( 7 )„ nV 

n>l 



and compare this with the sum 



E 



(SL)2 
V 2 In 



^4(n+l) ( 7 )„ nV 
We see that G and F are related by the expression 



(SL)2 
V 2 In 



G F + E 4n(n + 1) (7)n „ r 



(3.1) 



(3.2) 



(3.3) 



The new expression for the sum thus obtained will be easier to approximate since 
fewer terms will be required for a given accuracy. Moreover, using the Pochhammer 
identity (z) n +i = z(z + 1)„ , we note that F can be written in terms of a special form 
of the Gauss hypergeometric function [9] 

2 F 1 (a,b;c;z) = ^ 



n=0 



(c)„ n! 



(3.4) 



(with a circle of convergence |z| = 1 ) as 

(7-1) 



F = 



4(f -I) 2 



1; 7 -1;1 -1- 



(f-1) 2 
(7-1) 



(3.5) 



This generalizes the weak coupling expansion derived by Aguilera-Navarro and 
Guardiola to study Eq.(1.4) for a < 7 + 1 . However, we should note here the correct 
form of function F of the case 7 = 3/2 or A = that is 



F = 



1 



8(f -l) 2 



1 '2" 1; 2 ;1 



l-2(f-l) 2 



(3.6) 



not as quoted in [1-3]. Eq.(2.8) and (3.5) can be used now to obtain a sum for the 
infinite series 



E 

n>l 



f«)2 



4n(n+ 1) (7)„ n! 



— 4F 3 (1,1,- + 1,- + 1;7+1,2,2;1) 
I67 2 2 



(7-1) 



4(f-l) 2 



1; 7 -1;1 



(f-1) 5 



(7-1) 



(3.7) 
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valid for a/2. For the special limit a = 2 , the expression (3.5) has no meaning. 
However, the sum in the A 2 coefficient of (2.6) converges and follows from (2.8) by 
setting a = 2. Indeed, in this case, the sum in (2.6) becomes 



E 



(in 



= E 



{n-l)\ 



^ An ( 7 )„ n\ ^4 (7) 

n>l v ' n>l v 7 



E 

n=0 



(1) 



4 (7)71+1 n\ 47 



2 F 1 (1,1; 7 +1;1), 



(3.8) 



which follows immediately from (2.8). The Gauss hypergeometric function 2 -Fi(l, 1; 7+ 
1; 1) can be evaluated using the identity [9] 



2 Fi(a,b,c\ 1) 



to obtain 



E 



T(c)T(c-a-b) 
Tic-a^ic-h) 1 



c - a - 6 > 0, c> 6 > 



7 > 1. 



^4n( 7 )„n! 4( 7 - 1) ' 
Thus, for the case a = 2 , the weak coupling expansion (2.9) becomes 

E(a = 2) = 2VB, + J r) X-^ w X> + ..., 7 >1 

and finally, for B = 1 , A = or 7 = 3/2 , we get 

E(a = 2) = 3 + 2A - 2A 2 + . . . , 
as we expect. 

IV. More closed-form sums of infinite series 

For the special limiting case a — > 2 , we introduce a parameter e = f - 1 which will 
be chosen to approach zero. The function F , as given by (3.5), becomes in this limit 



(3-9) 
(3.10) 

(3.11) 
(3.12) 



(W - 11 
lim F = — lim e~ 2 

e^O 4 e^O 



2 Fi(e,e,7-l,l)-l 



(4.1) 



Using the series expansion of the Gauss hypergeometric function (3.4), Eq.(4.1) can 
be written, using T(z + 1) = zT(z) , as 



limF 



00 

iE 



r(n)r( 7 ) 



e^o 4 ^ nT(n + 7 - 1) 4 ' 



(4.2) 



Some results similar to Eq.(4.1) and (4.2), were first published, without detailed 
proofs, by Mitchell [11]. Now, using the identity (3.9), Eq. (4.1) can now be rewritten 

as 

T( 7 - l)r( 7 - 1 - 2e) 



lim F = ^ ^ lim er 2 

e^O 4 e^O 



- 1 



(4.3) 



r 2 ( 7 -l-e) 

Now employing Znojil's method [3], we can obtain closed form sums for other infinite 
series involving Gamma functions. Indeed, the Maclaurin's expansion of the gamma 
function 

r(c + x) = r(c)|l + xyj(c) + \[4> 2 + ^ (1) (c)] + • • • ] (4.4) 
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where ip(c) and V'"H C ) ,n>l, are the digamma and the polygamma functions [10], re- 
spectively. We can show, after expanding the gamma functions in (4.3) and employing 
multiplication and division of polynomials, that Eq.(4.4) can be written as 

lim^^VHl-l)^ (4-5) 
where ip( n \z) are the polygamma functions. Comparing (4.2) and (4.5), we have 

E^g^H^-). .>i. («) 

n=l y ' 

where tp (1 ^(z) is the trigamma function [7]. For the purpose principally of verification 
we now note some special cases. For 7 = 2 

00 1 2 

n— 1 

and for 7 - 1 = m > 2 positive integer, we find 

- T(n)T(m) £ _ " 1 

^nr(n + m) 6 ^ 2 (fc-l) 2 ' l " J 
by using the recurrence relation 

</,(")(> + 1) = V (n) (^) + {-i) n n\z- n -\ 

Further, for 7 = 3/2 

y r(n)r(i/2) ^ 

^nT(n+ 1/2) 2' 

and 

g r(n)r(mz) 1 g (1) fc 
nl h| mz w ! ' m 

n=l y ' k=0 

Finally, we can now have a finite sum for the infinite series (3.7) for the case a = 2 
and 7 > 1 , 

S 4n(,, + ( !)''h)„ 5 " £ »«(Mn + l; 1) - ^'"h " » + J- («» 



n> 



V. Conclusion 

We have obtained a compact weak-coupling expansion (2.9) for eigenvalues of 
the spiked harmonic oscillator Hamiltonian. Our expansion extends the earlier work 
of Aguilera-Navarro and Guardiola for 7 ^ 3/2, and it allows for arbitrary spatial 
dimension N and also, for N > 2, arbitrary orbital angular momentum I. Moreover, 
with the closed-form expressions we have been able to provide for the coefficient of 
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the A 2 term, the new expansion is easier to handle and calculate with, even at or 
near to the special value a = 2 . These analytic expressions describe approximately 
how the eigenvalues depend on all the parameters in the Hamiltonian. Such formulas 
are complementary to data obtained with the aid of a computer; moreover, they are 
useful in guiding a procedure that searches for very accurate numerical eigenvalues. 
As a byproduct of this work, we have been led to some simple closed-forms for a 
variety of interesting infinite series involving sums and ratios of gamma functions. 
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